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Evolution of a curved vortex filament into a vortex ring 


Parviz Moin, Anthony Leonard and John Kim 
NASA Ames Research Center, Moffett Field, California 94035 


The deformation of a hairpin-shaped vortex filament under self-induction and 
m the presence of shear is studied numerically using the Biot-Savart law It is 
shown that the tip region of an elongated hairpin vortex evolves into a vortex ring 
and that the presence of mean shear impedes the process. Evolution of a finite- 
thickness vortex sheet under self-induction is also investigated using the Navier- 
Stokes equations. The layer evolves into a hairpin vortex which in turn produces 
a vortex ring of high Reynolds stress content. These results indicate a mechanism 
for the generation of ring vortices in turbulent shear flows, and a link between the 
experimental and numerical observation of hairpin vortices and the observation of 
Falco [Phys. Fluids 20, S124 (1977)] of ring vortices in the outer regions of turbulent 
boundary layers 
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I. INTRODUCTION 


Hairpin-shaped vortices (also referred to as horseshoe-shaped vortices ) are com- 
monly proposed as the characteristic organized structures in turbulent boundary 
layers. 1-5 Recently, using vorticity fields generated by large-eddy simulation, Moin 
and Kim 6 demonstrated that turbulent channel flow contains an appreciable num- 
ber of hairpin vortices often inclined at 45° to the flow direction. It was also shown 
that hairpin vortices are associated with the bursting process, thus, indicating their 
importance in the turbulence energy production mechanism 7 . The same structures 
are found in homogeneous turbulent shear flow 8 and this implies that hairpin vor- 
tices are the characteristic structures not only in turbulent boundary layers, but in 
all turbulent shear flows. 

Falco 9 identified a relatively small-scale coherent motion, which he labeled “typi- 
cal eddy,” in the outer region of a smoke-filled, turbulent boundary layer . Using si- 
multaneous hot-wire anemometry and flow visualization, Falco 10 showed that these 
eddies produce most of the Reynolds stress in the outer region. Falco’s reported 
visualizations are generally in a plane parallel to the flow and perpendicular to the 
wall ( 2 : — y planes), but he does present sketches of the typical eddy in planes per- 
pendicular to the wall and the flow direction (y — z planes). Figure 1 shows Falco’s 
sketch of two views of the typical eddy. He points out that “sometimes only the 
solid outline (in Fig. la) is observed, but at other times the outline of the typical 
eddy consists of both the solid and dashed lines.” Because the cross-section outline 
(in an x — y plane) of a vortex ring inclined to the flow direction resembles the 
complete outline of the sketch in Fig 1(a), Falco identifies the typical eddy as a 
vortex ring It should be pointed out, however, that the cross-section outline of a 
hairpin vortex would resemble the solid line in Fig. l(a). 

A relevant computational work is that of Hama 11 who studied the deformation 
of an initially parabolic vortex filament by using the localized induction concept. 
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His calculations show that the lifted portion of the filament progressively assumes 
circular shapes resembling part of a ring vortex. Using the more accurate Biot- 
Savart induction law, Dhanak and De Bernardinis 12 demonstrated that an initially 
elliptic vortex ring of relatively high eccentricity breaks up into two smaller loops, 
but that an elliptic vortex of low eccentricity may simply oscillate without pinch-off. 

The common feature of the hairpin vortices observed in the simulations of Moin 
and Kim 6 was the bottle-neck, or U-shape, of their tip region. This observation and 
the results of Hama led the authors to speculate that the hairpin vortices observed in 
turbulent shear flows may pinch-off into ring vortices thus, providing a mechanism 
for the generation of ring vortices in turbulent shear flows and a link between the 
ring vortices observed by Falco and the hairpin vortices observed (or conjectured) 
by others. 

The objective of this study is to examine the motion of a curved vortex filament 
or a hairpin vortex under its self-induction and in the presence of mean shear. We 
were particularly interested in determining whether the hairpins pinch-off into ring 
vortices and, if so, what factors influence this process. This paper extends Hama’s 
work as follows: Hama used the local-induction approximation, which neglects the 
contribution to the velocity from other portions of the vortex. In the present work 
the velocity is calculated using the full Biot-Savart line integral and also by means of 
three-dimensional Navier-stokes equations. The Navier-Stokes computations were 
performed, in part, to account for viscous effects which are significant when vortices 
of opposite sign merge in the pinching process. Hama’s work did show the evolution 
of a curved vortex filament into the fi-shape, but his calculations were not carried far 
enough to indicate the possible occurrence of pinch-off. The filament may oscillate 
about a configuration as has been shown for an elliptic vortex ring by Dhanak and 
De Bernardinis . 12 In any case, it is doubtful that the local- induction approximation 
could lead to meaningful results when the vortex cores approach each other in the 
splitting process. One may consider three stages in the pinching process . the early 
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stage where the local induction effects are dominant, the stage during which the 
influence of the neighboring vortices is important, and the final period where viscous 
stresses become significant. The Biot-Savart calculations encompass the first two, 
and the Navier-Stokes calculations encompass all three. We will also present the 
effects of a background shear flow and the presence of a wall on the evolution of the 
vortex filament. Finally, we will use our results to propose a model for the structure 
of the vorticity field in turbulent shear flows. 


II. NUMERICAL METHODS 

In the Lagrangian (Biot-Savart) calculations, we assume that the vorticity field 
is in the form of a tube of concentrated vorticity following the space-curve r(£, t) 
plus weak background vorticity producing a mean shear. The motion of the space 
curve is determined by self-induction plus an external component, 


dr _ _T_ f 

dt 4tt J 


dr _ r f [r({,t) ~r(C,i)l x jfig(|r-r'[) 


w ]3 d(' + V(r,t) 


r — r 


( 1 ) 


where U is the background velocity, T is the circulation of the vortex tube, and the 
function g models the effect of distributed vorticity within a core of effective radius 
a (see Ref. 13). The equation of motion is valid in the thin-filament approxima- 
tion: the structure of the core remains nearly constant in time and the disturbance 
wavelengths are larger than the core radius. In addition, we neglect the generation 
of new disturbance vorticity produced by the velocity field of the hairpin vortex 
acting upon the background vorticity (see Sec. IIIC). 

For numerical purposes, the space-curves were tracked by following a sequence of 
node points on the curve. The function g was taken to be 


g{y) = 




(y 2 + acr 2 ) 3 / 2 
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with a = 0.413 (Ref. 14). Integration along the curve was performed by the 
midpoint rule and second order central differencing was used to estimate J-fr. The 
time-integration method was third-order Runge-Kutta. 

The spacing between the node points, A, was sufficiently small so that length 
scales of interest would be represented accurately. In particular, if the initial mini- 
mum radius of curvature of the hairpin vortex is p 0 , then 

A < p 0 (2a) 

On the other hand, the cutoff scheme used in the thin-filament approximation gen- 
erally produces spurious results for length scales smaller than a. The constraint 

A > 2 a (26) 

was also enforced to prevent the appearance of these scales. 14 

A pseudo-spectral method was used for the Navier-Stokes computations. 15 The 
flow is assumed to be between two parallel plates and is maintained by a mean 
pressure gradient. No-slip boundary conditions are enforced on the plates, and 
the flow is periodic in the streamwise and spanwise directions. In the direction 
normal to the walls, y, the flow field is represented by a series of 129 Chebyshev 
polynomials. Fourier series with 128 modes are used in each of the streamwise and 
spanwise directions. Thus, a total of about 2 x 10 6 collocation points were used. 
Time is advanced by a second-order semi-implicit scheme. 


III. RESULTS AND DISCUSSION 

The results presented in sections IIIA-IIIC are obtained using the Biot-Savart law 
and those in Sec. HID are obtained using the Navier-Stokes computations. The flow 
is in the x-direction, y is the direction normal to the wall(s), and z is the spanwise 
direction. 


5 



A. Evolution of an isolated parabolic vortex 


The initial configuration of the vortex filament is the parabola x = — az 2 ; y — 
0. The constant a has the dimension of inverse length. This filament is isolated 
in space; that is, the effects of shear or image vortices accounting for the wall 
impermeablity are not included. Note that in our calculations [when U — 0 in Eq. 
(l)] the actual value of the circulation T is not significant, since it can be absorbed 
into a re-scaling of the time-variable; the only nondimensional length parameter is 
era. We are assuming that our vortex filament is formed from the roll-up of a sheet 
of vorticity in a turbulent boundary layer (see section HID). The diameter of the 
vortex core is specified to be <7 = 0.125. In Fig. 2, the evolution of the filament 
with o = 1 is shown. At each instant in time, the plan, side and end views are 
presented. The range of the variable z was from z = —20 to z = 20, with x ranging 
from 0 to -400. The large range of the independent variables is used to minimize 
the effects of artificially terminating the vortex filament on the deformation of the 
tip region. It will be shown that such deformations are confined to Az ~ 50 in the 
tip region and hence unaffected by the end conditions. The filament is represented 
by N = 129 elements uniformly spaced in the ^-direction. The relatively small core 
size and the value of N are chosen in accordance with the accuracy requirements 
[Eqs. (2a), (2b)]. 

The small radius of curvature in the tip region (p = ^ = 0.5) results in the rapid 
rise of the tip above the x — z plane, giving a three-dimensional configuration to the 
vortex. Gradually, other portions of the vortex in the vicinity of the tip also rise 
It appears that the raised portion has approximately uniform velocity in the x — y 
plane. This is because the subsequent cross sections of the filament in this plane 
are approximately parallel to each other. As the tip region lifts from the plane 
of the initial filament, secondary curvature is introduced between the legs and the 
tip. This curvature in x — y plane induces a z-velocity component which forces the 
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legs away from each other and leads to the bulged-out shape, or f2 shape, evident 
in the plan and end views of Fig. 2. The curvature in the x — y plane and the 
subsequent velocity in the 2 -direction create a curved section in the x — 2 plane 
leading to negative velocity in the y-direction and to the corresponding downward 
movement of the filament in the transition regions between the legs and the tip. 
The resulting curvature in the x — y plane induces opposite velocities on each leg 
in the 2 -direction, forcing them toward each other and causing pinch-off. The plan 
and end views show that the tip region progressively evolves into a closed loop 
resembling a ring vortex. 

The calculations are terminated at the nondimensional time ta 2 T = 6.5 x 10 4 
when pinch off is deemed imminent, and the spatial accuracy of the computations 
has begun to deteriorate. When the cores of vortices of opposite sign approach each 
other, viscous diffusion and other fine scale effects become significant. Our inviscid 
calculations have no provisions for accommodating viscous effects. In the absence 
of viscous effects and during the final stages of the pinching process, high, localized 
induced velocities lead to increased vortex stretching and increased local curvature. 
In fact, numerical studies by Siggia 16 suggest that the final stages of the pinching 
process at high Reynolds numbers involve intense regions of vortex stretching in 
which the self-induced motion of opposite-signed vortex filaments generates rapidly 
increasing filament line lengths, requiring the insertion of a large number of addi- 
tional computational node points. We see what might be the beginning of such a 
cascade process in the above Biot-Savart calculations, so that there are strong indi- 
cations that pinching does occur. In this investigation, Navier-Stokes computations 
(see Sec. HID ) were performed in part to enable simulation of the final stages of 
the pinching process. 

In what follows, we will refer to the time required for the two legs of the filament 
to most closely approach reach other (as allowed by the accuracy of the computa- 
tions ) as the pinching time. Computations with larger core sizes (not shown) could 
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be extended to a point where the two cores actually overlapped without encounter- 
ing high local curvature in the filament. Note that our estimate of pinch-off time is 
admittedly crude. Because of numerical errors encountered during the pinching pro- 
cess, this estimate is probably dependent on numerical method and on the number 
of nodes used. However, in this study we are mainly interested in whether pinching 
occurs. For further examination of the process itself, more refined techniques, such 
as that used by Siggia 16 should be resorted to. 

For isolated parabolic filaments, the only nondimens ional length scale present 
is era. Experiments with different values of era show that except for a time-scale 
dependence, the pinching process takes place independently of oa and that the 
nondimensional time for pinch-off depends weakly on oa. Note that the velocity 
of the filament owing only to local induction effects is proportional to ln(^). The 
evolution of the parabolic filament with a = 0.4 and a = 0.0625 (oa = 0.025) is 
shown in Fig. 3. For this calculation, N = 321 node points were used to better 
resolve the curvature of the filament in the tip region while maintaining the accuracy 
requirement [Eqn. (2b)]. To achieve geometric similarity with the previous case, 
the range of the variable z was from z = —50 to z — 50. Once again the filament 
appears to split into a ring vortex (at ta 2 T = 3.4 x 10 4 ). It is interesting that 
due to the high curvature of the filament near the pinching region, a second ring 
begins to form just before the section where the first ring is pinched away. The 
same tendency might have occurred in Fig. 2 had more nodes were used. In this 
calculation, approximately 8.6 x 10 5 time-steps were required. The slight asymmetry 
developed in the shape of the filament is due to the accumulation of round-off errors. 

The above results indicate that the tip region of an isolated parabolic vortex 
filament resembling a stretched hairpin vortex evolves into a ring vortex. 
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B. Effect of image vortices 


To account for the wall impermeability, an image vortex is introduced in to the 
calculations which will also introduce another length scale. We will consider the 
same parabola as that in Fig. 2 (x = — z 2 , a = 0.125) placed above the wall 
at y = 5. The image vortex is at y = —5, and each vortex contains N = 129 
elements. The evolution of this filament is shown in Fig. 4. The main effect of the 
impermeable wall or of the image vortex appears to be enhancement of the pinching 
process. The velocity induced by the image vortex moves the legs of the parabola 
toward each other. This effect increases as the filament is moved closer to the 
wall and its image. Note that similar reasonings would indicate that a commonly 
assumed configuration of the wall layer of turbulent boundary layers consisting 
of elongated pairs of counter-rotating longitudinal vortices is essentially unstable. 
That is (depending on whether the flow between the vortices is toward or away from 
the wall), counter-rotating longitudinal vortices near a wall either merge or move 
apart. 


C. Effect of mean shear 

We approximate the effect of mean shear on the vortex filament by adding the 
mean velocity profile U{y ) to the right-hand side of the 2 -component of Eq. (1). 
As pointed out by Hama 17 and discussed by Aref and Flinchem 18 , this treatment 
of the background or mean shear is an approximation which neglects the effect of 
the filament on the background flow. In particular, we neglect the effect of new 
perturbation vorticity which is estimated to be generated at a rate proportional to 
(u/L)fi, where H = |^| is the background shear, u is the characteristic (large-scale) 
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velocity induced by the hairpin vortex, and L is the large length scale of the vortex. 
By the time pinch-off occurs, t = T, the neglected vorticity has the magnitude 

Acj « (y)nr 

1j 

But ( Vf/L)T — 0(1), where Vf = 0[ufn(£)] is the characteristic velocity of the 
hairpin vortex filament. Therefore, 

Au 1 

Ml) 

and as long as the core radius is much smaller than the characteristic radius of 
curvature, the approximation appears valid. 

We shall first consider the case of uniform shear U(y) — 0.2 y and the isolated 
filament (x = — z 2 ,a — 0.125). The results are presented in Fig. 5. When Fig. 
5 is compared with Fig. 2, it is clear that the effect of uniform background shear 
is to inhibit or retard the pinching process. This can be explained by noting that 
the points with higher y-locations on the filament are moved at a faster rate in the 
x-direction than those at lower y- values, thus restraining the tip region from taking 
a circular shape. The effect of stretching by the mean flow is evident; the node 
points in the tip region are rapidly moved apart. 

Next, instead of uniform shear, we impose a boundary-layer-type profile, 

U(y) = c[ 1 - (0.02y - l) 6 ] y < 50 

U(y) = c y > 50 

on the filament. With c = 3.75, this profile has the same shear rate at y = 7.5 
(where the filament is placed initially) as in the uniform shear case; but it has 
smaller shear rates at higher y-locations. The evolution of the filament is shown 
in Fig. 6. The lifted portion of the filament is exposed to lower shear rates, and 
hence the restraining effect of shear is less pronounced. In agreement with Falco’s 
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observations, these results indicate that it is more likely that vortex rings would be 
detected in the outer region of turbulent boundary layers. It is emphasized that 
the presence of mean shear is essential for the formation of elongated vortex loops 
or hairpins with which the calculations were initialized; but the resulting hairpins 
more readily pinch-off into vortex rings in regions with small mean shear. 

D. Evolution of a vortex sheet 

The results presented so far were obtained from inviscid Biot-Savart integral 
calculations which did not account for the effect of viscosity. The actual merging of 
the two legs of the hairpin requires a calculation that includes viscous effects. It is 
also desirable to initialize the computations with a flow structure that occurs often in 
turbulent boundary layers. In shear flows, the undisturbed mean flow is a continuum 
of vorticity which is in the form of an infinite sheet of finite thickness. It has been 
conjectured that a perturbation of this sheet by random velocity fluctuations causes 
it to roll up into a vortex loop which when stretched appears as a hairpin vortex. 6 

N 

To follow the evolution of an isolated hairpin vortex from its inception, a three- 
dimensional, time-dependent calculation of the Navier-Stokes equations was carried 
out. The flow geometry is a channel flow with a constant mean pressure gradient. 
The Reynolds number based on the centerline mean velocity and channel half- width 
is about 3500. The initial condition consists of the mean-velocity profile of a fully 
developed channel flow plus a perturbation in the form of a deformed vortex layer 
in the lower half of the channel. Using the shear velocity u r , obtained from this 
mean-velocity profile, the wall-coordinate of the channel centerline is calculated 
to be = 6u T lv = 180. The computational grid resolutions in the streamwise 
and spanwise directions are Ai + = 18 and A z + = 6, respectively. In the normal 
direction, nonuniform mesh spacing is used which ranges from 0.1 wall units near 
the wall to 4 wall units near the channel centerline. 
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The vortex lines drawn in the vicinity of the initial perturbation are shown in 
figure (7a). Note that the mean velocity field contains only the spanwise vorticity 
component, and that departure of the vortex lines from ani-z plane represents a 
perturbation. Figures 7(b)-7(d) show the evolution of this perturbed vortex layer. 
In contrast to the Biot-Savart calculations reported above, vortex lines are traced 
using the computed velocity field until they intersect the boundaries of the domain 
or are trapped in a loop 19 It can be seen that near the vortical structure, vortex 
lines are stretched and coalesce into vortex tubes, indicating that the initial vortex 
sheet has rolled up into a hairpin vortex. Figure 7(c) shows that the induced motion 
of the hairpin vortex draws the downstream vortex lines toward itself. Recall that 
local curvature of the filaments induces velocities in the z-direction that produce 
the characteristic bulged-out, or Cl, shape. When the downstream vortex lines 
are drawn toward the hairpin vortex, these velocities are enhanced because of the 
superposition effect, thereby accelerating the pinching process. As shown in Fig. 
7(d), these filaments eventually pinch-off to form a vortex ring. 

The formation of the vortex ring just described is slightly different from the one 
we discussed in the previous sections: in the previous examples the hairpin vortices 
themselves deformed into vortex rings, whereas in this example the hairpin vortex 
creates a vortex ring from the surrounding vortex lines. In both cases, however, the 
rings are formed because of the induced motion of a curved vortex filament or a 
hairpin, as discussed in sec. IIIA. Examination of the contour plots of the velocity 
product uv indicates that the vortex ring shown in Fig. 7(d) is associated with high 
Reynolds shear stress, which is consistent with the measurements of Falco. 10 The 
point of maximum — uv is located just downstream of the ring. 


12 



IV. CONCLUSIONS 


We have demonstrated by numerical experiments that a curved filament of con- 
centrated vorticity evolves into a vortex ring as a result of self-induction effects. 
This result suggests a mechanism for the generation of ring vortices in turbulent 
shear flows. The ring vortices are more readily generated in regions with relatively 
low mean shear. A three-dimensional, unsteady calculation of the Navier-Stokes 
equations showed that a finite-thickness layer of vorticity, resembling a perturba- 
tion to the mean flow in a turbulent boundary layer, rolls up into a hairpin vortex 
which in turn produces a vortex ring having high Reynolds shear stress. These 
results provide a link between Falco’s observation of ring vortices in the outer re- 
gion of turbulent boundary layers and the presence of hairpin vortices reported by 
others. 

The observations of hairpin and ring vortices and the results of the present in- 
vestigation suggest that in turbulent boundary layers, perturbation of a layer of 
vorticity leads to its roll-up into filaments of concentrated vorticity. These fila- 
ments are then stretched into hairpin vortices. Depending on their lifetimes, the 
tip region of some of these hairpins may evolve into vortex rings. 
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A selection of these lines is shown in Figs. 7(c) and 7(d) In Fig. 7(c) none 
of these lines formed a loop, whereas in Fig. 7(d) several lines initiated from 
different grid points formed vortex loops. 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Sketches of the organized structures in the outer layer by Falco 10 (a) Side view, 
(b) front view. 

Evolution of the parabolic vortex filament x = — z 2 , o = 0.125. Plan views are 
in the left column, the middle column shows the side views and the end views 
are shown in the right column: ia 2 r = (a) 0; (b) 7.1 x 10 3 ; (c) 2.8 x 10 4 ; (d) 

5.0 x 10 4 ; (e) 6.5 x 10 4 . The symbols indicate the location of the node points. 
Evolution of the vortex filament x = — 0.42 2 , a = 0.0625. Plan views are in 
the left column, the middle column shows the side views and the end views are 
shown in the right column: ta 2 T = (a) 1.9 x 10 3 ; (b) 1.2 x 10 4 ; (c) 2.7 x 10 4 ; (d) 
3.4 x 10 4 . The symbols indicate the location of the node points. 

Evolution of the filament in Fig. 2 initially placed at y = 5 and the image vortex 
placed at y = —5: ta 2 T = (a) 7.1 x 10 3 ; (b) 2.1 x 10 4 . 

Evolution of the filament in Fig. 2 in the presence of uniform shear. ta 2 T = (a) 

7.1 x 10 3 ; (b) 2.1 x 10 4 ; (c) 2.8 x 10 4 . 

Evolution of the filament in Fig. 2 in the presence of a boundary-layer- type 
background shear. ta 2 V — (a) 2.8 x 10 4 ; (b) 5.0 x 10 4 ; (c) 7.5 x 10 4 ; (d) 7 8 x 10 4 
Evolution of a perturbed vortex layer. tu T jS = (a) 0; (b) 1.2; (c) 1.56; (d) 2.16. 
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